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ABSTRACT 

Context. On March 2010, the symbiotic binary V407 Cyg erupted as a result of a nova explosion. The event gave rise to a two- week 
long burst of >100MeV y-rays detected by FermifLKI, a unique observation testifying to particle acceleration in the system. 
Aims. The outburst can be considered a scaled-down supernova, with short dynamical time scale, and thus can constitute a test case 
for theories of the origin of galactic cosmic rays. We aim at determining the properties of the accelerated particles and identifying the 
origin of the high-energy radiation. 

Methods. We developed a model for diffusive shock acceleration and non-thermal emission in V407 Cyg, complemented by an eval- 
uation of the thermal emission from the shocked plasma. We considered both leptonic and hadronic contributions to the non-thermal 
processes, and investigated the effect of many binary and nova parameters. 

Results. The y-ray emission is mostly of leptonic origin and arises predominantly from inverse-Compton scattering of the nova light. 
Upscattering of the red giant photons is a minor contribution. Matching the light curve requires gas accumulation in the vicinity of 
the white dwarf, as a consequence of wind accretion, while the spectrum imposes particle scattering close to the Bohm limit in the 
upstream equipartition magnetic field. The nova accelerated protons (respectively electrons) with energies up to ^300 GeV (respec- 
tively =:20GeV), for a total non-thermal energy =d0 43 erg after two weeks, representing =T0% of the initial nova kinetic energy. The 
electron-to-proton ratio at injection is 6%. 

Conclusions. The V407 Cyg eruption can be understood from the same principles that are invoked for particle acceleration in super- 
nova remnants, although without the need for strong magnetic field amplification. The population of novae in symbiotic systems is a 
negligible source of Galactic cosmic rays, and most likely not a class of TeV-emitters. Together with classical novae, their high-energy 
emission provides a new perspective on particle acceleration at shocks. 

Key words, binaries: symbiotic - novae, cataclysmic variables - stars: individual: V407 Cyg - acceleration of particles 
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1. Introduction 

V407 Cyg is a binary system belonging to the class of symbiotic 
binaries, in which a hot compact component is accreting matter 
from a late-type giant predominantly via a stellar wind. These 
systems are characterised by composite optical spectra show- 
ing absorption feature s typical of late-type g iants and strong 
nebular emission lines dBelczvnski et al.lt2000l) . V407 Cyg con- 
sists of a white dwarf (WD) and a Mira-type M6 III red giant 
(RG) with a 745-day pulsation period, which places the sys- 
tem in the group of D-type (dusty) symbiotic stars or symbiotic 
Miras (Munari et al.lll990T) . From the RG pulsation period and 
Mira period-luminosity relation, the distance to V407 Cyg is es- 
timated at 2.7 kpc. Isolated Mira variables with such long pul- 
sation periods generally have thick circumstellar dust envelopes 
detectable in infrared but causing high obscuration in the opti- 
cal. In V407 Cyg (and more generally in symbiotic Miras), how- 
ever, the influence of the WD through its orbital motion, ener- 
getic radiation, and repeated outbursts inhibits dust formation 
over most of the Mira wind except in a shadow cone behind the 
RG. As a consequence, V407 Cyg can be observed in optical, 
with characteristic obscuration events due to the dust cone. An 
orbital period of about 40 yr was inferred from the dust obscu- 
ration events, which for a RG mass of 1.0 M Q and a WD mass 
of about 1.2 M p correspon ds to an orbital separation of about 
16 AU dMunari etalJfl990h . 

V407 Cyg has been monitored in optical for about 75 years, 



although irregularly. Apart from the long-term modulation of 
the Mira pulsed emission by dust obscuration, two episodes of 
enhanced activity were observed in 1936 and 1998 (the first 
one being the eve nt that led to the discovery of V407 Cyg; see 
lHoffmeisterf [T949). On these occasions, the optical emission of 
the WD rivaled that of the Mira star at its maximum and then 
slowly declined over the following several years. These flares 
are typical of all symbiotic bi naries and probably r esult from 
changes in the accretion rate dKolotilov et alJ 120 03). In 2010, 
V407 Cyg was observed on March 10.813 (UT) at V=7, about 
5 magnitudes above the maximum of the Mira cycle, a bright- 
ness level never reached before in the photometric history of 
the system (Nishivama et aljfeOlOl) . Early specttoscopic analy- 
ses described the event as a He/N nova expanding in the wind 
of the Mira com panion, with some s imilarity with the 2006 out- 
burst of RS Oph (M unari et al.ll2.Ql6T) . Follow-up observations of 
V407 Cyg were performed in various bands including radio, in- 
frared, X-rays, and y-rays. In particular, the nova explosion was 
observed in y-rays >100MeV by FermifL AT, and it was the first 
detection of such an object at high energies. 

In the present paper, we focus on the non-thermal emission 
of V407 Cyg, especially that detected by FermifLAT. This radia- 
tion testifies to rapid particle acceleration in the nova blast wave 
as it propagates into the RG surroundings. Whereas in classical 
novae shock and ejecta are thought to propagate almost freely in 
a tenuous ambient medium, in symbiotic systems they interact 
with a denser environment and are thus slowed-down more effi- 
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ciently as larger quantities of material are swept up and heated 
up (thermally or non-thermally). This makes V407 Cyg a short- 
lived particle accelerator, which in some respect behaves like a 
scaled-down supernova remnant with a life time of just a few 
weeks or months. The system therefore can constitute a valu- 
able test case for theories of the origin of galactic cosmic rays. 
We present a model for diffusive shock acceleration and non- 
thermal emission in V407 Cyg, complemented by an evaluation 
of the thermal emission from the shocked plasma. We consid- 
ered both leptonic and hadronic contributions to the non-thermal 
processes, and investigated the effect of many binary and nova 
parameters. From a comparison to the various observables, we 
derived estimates for the maximum particle energies, the non- 
thermal energy budget, the relative contribution of electrons and 
protons to the emission, and the parameters of the binary system 
and nova event. 

2. Observations 

In contrast with its previous flares of 1936 and 1998, the outburst 
experienced by V407 Cyg in March 2010 was a thermonuclear 
runaway, the violent burning of a shell of accreted material in 
degenerate conditions. The high temperatures reached cause the 
rapid expansion of the shell and eventually its partial ejection 
from the WD, at velocities of order 1000 km s~'. Emission line 
widths provide information about flow velocities after the out- 
burst. Ha observations on 13 March revealed a line width consis- 
tent with an expansion velocity of the emitting material of 2760 
kms^FWHM, full width at half maximum). Such a broad line 
originates either from the material ejected in the nova eruption 
or from the post-shock medium (through charge-exchange), and 
the inferred velocity at day 2.3 after optical detection is a lower 
limit but most likely a close estimate of the initial speed of the 
ejecta. Subsequent line observations revealed that the bulk of the 
emitting material is progressively slowed down, with an FWHM 
expansion velocity droppi ng to 1500km s~' a t day 6.3 and then 
to 400 km s 1 at day 48.2 dMunari et alJl201 lh . Yet, the extreme 
radial velocities of the Ha line profile are still about 2500 km s -1 
even at day 48.3 after optical dete ction, suggesting th at parts of 
the ejecta were hardly decelerated (Shore et al. 2 01 lh . 

The outburst in >100MeV y-rays started on the second half 
of 10 March 2010, the same day as that of the optical maxi- 
mum (note that the exact start of the nova is uncertain by up to 
3 days because the last pre-outburst optical observation was car- 
ried out on 8 March). The flare lasted for two weeks and a peak 
was reached between 13 and 14 March with a flux a factor of 
2 higher than that of the initial detection. The average spectral 
energy distribution over the active period is a power law with a 
significant exponential cutoff, and a total flux >100MeV of (4.4 
+0.4) x 10 _7 phcm _2 s _1 . No evidence for spectral variability 
was found over the active y-ray period. 

In the radio band, measurements at 30 GHz were made by the 
OCRA collaboratiorQ on several epochs between 24 March 2010 
and 3 October 2010 . They revealed an emission at the 25 mJy 
level around day 15, rising to a peak flux density of about 45 mJy 
at day 40, and then followed by a drop back to about 25 mJy at 
day 100. Other observations came from the JVLA Nova TearrQ, 
at several frequencies between 1 .4 and 45 Ghz and covering a 
period from day 24 to 470. These data show that the radio max- 
imum occurs later and lasts longer at lower frequencies: at 1- 
2 Ghz, a flat maximum at 6-7 mJy is observed over day 50 to 

1 See http://www.jb.man.ac.uk/research/ocra/v407cyg 

2 See https://safe.nrao.edU/evla/nova/#v407cyg 
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Fig. 1. Illustration of the density structure adopted for the 
medium through which the shock travels. It comprises a 
spherically-symmetric wind profile about the RG and possibly 
a Gaussian-shaped matter accumulation around the WD. 



150 approximately. The spectrum of the radio emission flattens 
progressively over time, from about F(v) oc y - 7 at day 24 to an 
almost flat spectrum at day 360. 

X-ray emission from V407 Cyg was detected with 
Swift/XKT in the 0.3-10keV band. Observations started 3 days 
after the optical maximum and lasted for about 4 months with 
recurrent visits every few days. Over that period, the emission 
rose rapidly from about day 10, reached a maximum at day 30 
with a luminosity Lx — 1 - 2 x 10 34 erg s _1 , declined slowly un- 
til day 60, and then remained appro ximately constant at about 
half the peak luminosity until day 90 dShore et alj|201 lh . The X- 
ray spectrum consists of two components: a soft blackbody and 
hard, optically thin emission from a high temperature plasma 
dNelson et alj l2012h . The blackbody component is consistent 
with the emergence of supersoft emission from the nova whereas 
the hard X- rays trace shock-heate d material, as previously seen 
in RS Oph (So koloski et alj|2006h . 

In the following section, we describe our model for the 
nova outburst and accompanying non-thermal processes. The 
key facts that we sought to reproduce are the early appearance 
and rapid drop of y-ray emission over days 0-14, the level of 
the thermal emission over days 20-40, and the effective slowing- 
down of at least a fraction of the shock front at day 40 after 
outburst. As demonstrated below, thermal and non-thermal con- 
straints should be considered simultaneously because they do not 
necessarily lead to the same requirements. 

3. The model 

3.1. Hydrodynamics 

The nova shell ejection drives a shock wave in the medium, 
where the sound speed is about a fewkms -1 . The discontinuity 
associated with the shock front is a fundamental ingredient of 
the first-order Fermi acceleration process, and our modelling of 
the nova hydrodynamics therefore focused on the shock position 
and velocity. The nova shell and associated blast wave were as- 
sumed to follow the typical behavio ur of a simple non-radiativ e 
supernova remnant (see for instance [Truelove & McKee1 [l999). 
In such a case, the early evolution is usually divided into two 
stages: first an ejecta-dominated stage (hereafter ED-stage), 
where the ejecta expand freely and the shock moves at a constant 
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speed, both being largely unaffected by the surrounding medium; 
second a Sedov-Taylor stage (hereafter ST-stage), where the ma- 
jority of the ejecta kinetic energy has been transferred to the 
swept-up ambient gas and the flow evolves adiabatically. For 
each phase, we adopted the limiting self-similar solutions and 
neglected the progressive non-self-similar transition in between. 
Instead, the change from ED-stage to ST-stage was set to happen 
when the swept-up mass equals the ejecta mass. 

While in classical novae the shell essentially expands freely 
into a tenuous circumstellar medium, in V407 Cyg the ejecta 
have to expand into a dense environment owing to the thick red 
giant wind. The nova shell can therefore be efficiently slowed 
down, especially the parts moving in the direction of the RG. 
The blast wave in V407 Cyg rapidly becomes aspherical because 
of the anisotropic matter distribution at the WD position. We 
therefore adopted a cylindrical geometry with origin at the WD 
position and symmetry in azimuth, about the WD-RG axis. Any 
point is characterised by an angle 9 from the symmetry axis and a 
distance r from the WD (with 9 = 0° being the direction from the 
WD to the RG). In our most basic model, the WD is embedded 
in the RG wind, which is assumed to be spherically-symmetric 
and to have constant velocity out to large distances. Yet, because 
of the accretion and orbital motion of the WD, the matter dis- 
tribution before outburst may be more complex than the simple 
wind profile, with a ccumulation of gas in the orbital plane and 
around the WD (see lWalder et alj|2008l for the case of RS Oph). 
We therefore also considered a possible circumstellar density en- 
hancement (CDE), although still in the 2D approximation. In the 
most general case, the density p tota i in the surrounding medium 
consists of two terms: 



p H ,(9, r) = p wf 



1 



r 2 + d 2 , — 1rd or u cos 



Pcde(Q, r) = pcde,o exp 



orb ' 

r sin 6 

bcDE 



Y COS 9 



CDE 



(1) 



(2) 



where p w (9, r) and pcde(&, r) are the wind and CDE density, re- 
spectively. The wind profile is defined by its normalisation, p w o, 
which is obtained from the stellar mass-loss rate M and constant 
wind velocity V w . The CDE profile is defined by its normali- 
sation, pcde,q, and by its characteristic scale lengths, Icde and 
bcDE- The parameter d or b is the orbital separation. The form of 
the CDE was taken from a recent study of the X-ray thermal 
emission from V407 Cyg (see Sect. 16. lb . Modelling of the nova 
optical lig ht curve also s uggests the presence of a large equato- 
rial disc dHachisu & Katoll2012l) . 

The shock wave evolution was computed from the following 
expressions for the shock velocity V s in direction 9 and at time t: 



V s (6,t) = V ej 

v s (e,t)= J 



for M, u < 



4tt 



M ei V 2 . 

ej ej 



for M,u > 



M, 



4-n 



(3) 
(4) 



where M e j and V e j are the mass and initial velocity of the ejecta. 
M su is the swept-up mass per unit solid angle, a quantity that 
increases with time but at a rate depending on the direction: 



M su (9,f) = 



2n sin# 
with R s {9, t) 



i r R ° 

n6d9 J Rw[ 

f 

Jo 



2np lotai (9,r)r 2 sm9&9&r (5) 
V s (9,t)dt (6) 



where R s is the shock radius and Rwd the WD radius. An implicit 
approximation in the above equations is that the shock structure 
has a negligible thickness and can be characterised at any (9, f) 
by a single position and a single velocity. 

3.2. Particle acceleration 



Our m odel for particle acceleration is inspired from lBall & Kirkl 
(fl99l and includes two zones where non-thermal particles are 
confined: an acceleration zone (zone A) fed by a constant frac- 
tion of the freshly swept-up material and where particles un- 
dergo diffusive shock acceleration (DSA) and competing energy 
losses; a cooling zone (zone B) fed by the downstream advection 
of particles from the acceleration zone and where particles un- 
dergo energy losses only. We use a thin-shell approximation and 
like the hydrodynamic shock structure, both zones are assumed 
to be coincident with the forward shock. This is a simplification 
especially for the cooling zone, since advected particles can in 
principle spread over a certain length in the downstream flow (al- 
thou gh that length rem ains a small fraction of the shock radius, 
see ISchure et alJuOlOh . In our model, they remain trapped just 
behind the forward shock and experience the immediate post- 
shock physical conditions. In addition, acceleration is assumed 
to occur at the forward shock only (no contribution from the 
reverse shock). The calculations are done in the test particle ap- 
proximation, which means that we neglected the non-linear feed- 
back effects of the non-thermal population on the shock structure 
and acceleration process. 

In the acceleration zone, charged particles with sufficient ki- 
netic energy can cross the shock multiple times and are scat- 
tered by magnetic perturbations on either side of the forward 
shock. Because of the velocity discontinuity, particles experi- 
ence a Lorentz boost at each shock crossing before being almost 
isotropised in the local fluid frame. The average rate of momen- 
tum gain for the diffusive shock acceleration process is 



dp\ 

dt I DSA 



V: 



g-1 

3q ) D(l +q)- 



(7) 



where q is the shock compression ratio (which is computed at 
all (9, f), but remains equal to 4 over the majority of the shock 
surface), p is the particle momentum, and D is the spatial dif- 
fusion coefficient. In the most general case, the diffusion coef- 
ficient has a dependence on space and momentum coordinates. 
In our model, we assumed that diffusion properties do not vary 
over the extent of the accelerator, hence the use of a single coeffi- 
cient D that applies both upstream and downstream. This implic- 
itly means similar magnetic field conditions on either side of the 
shock. In particular, the field is assumed to be highly tangled ev- 
erywhere (because of strong preexising turbulence upstream), so 
that we can neglect the issue of parallel versus perpendicular dif- 
fusion as a function of the propagation angle (which is thought 
to be important in some objects, e.g. SN1006; see IVolk et al.l 
120031 

The diffusion coefficient is one of the critical parameters of 
DSA and largely controls the maximum particle energy that can 
be attained. Particle scattering must be efficient enough to en- 
sure a large momentum gain rate. That it is the case in many 
supernova remnants is now supported by both observational ev- 
idence and theoretical developments. Magnetic field amplifica- 
tion by cosmic-ray streaming upstream should provide sufficient 
magnetic turbu lence at the rel evant scales for efficient resonant 
scattering (see ICapriolil |20 1 ll and references therein). Particle 
diffusion is then expected to occur close to the Bohm limit in the 
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amplified field, where the pitch angle scattering mean free path 
is of order of the gyroradius. In our model, we used a diffusion 
coefficient that is a multiple of the Bohm diffusion coefficient in 
the upstream non-amplified field: 



D (p) = C^nr = C D Bohm 

3eB s 



(8) 



where e is the elementary charge, B s the upstream non-amplified 
magnetic field strength just ahead of the shock front, and B — v/c 
where c is the speed of light. The factor g, which will be termed 
diffusion efficiency in the following, actually hides all the de- 
tails of magnetohydrodynamic wave amplification, advection, 
and damping. In the most general case, it is time-variable and 
evolves as the shock speed and upstream medium conditions 
change. For simplicity, we considered that ( is constant and 
treated it as a free parameter. Overall, with the adopted hypothe- 
ses, the momentum gain rate is independent of p. 

DSA is accessible to particles that have a large enough gy- 
roradius to probe both sides of the shock. In the thermal leakage 
scheme, only a small fraction ~ 10~ 2 - 10 4 of the downstream 
thermal p ool can enter the acceleration process as suprathermal 
particles (IVolk et al.ll2003l) . In that case, electrons are almost 
excluded because of their lower mass, hence lower momentum 
even for instant post-shock temperature equilibration; far less 
electrons are thus expected to be accelerated, and those that are 
have most likely been preacceler ated to relativis tic energies by 
some microphysical mechanism (Morlino 2009). Estimates of 
the electron-to-proton ratio from studies of non-thermal emis- 
sion from supe rnova remnants are in the range ~ 10~ 2 - 10~ 5 
dMorrinoll2009l) . while it is ~ 10~ 2 from measurements of the 
cosmic -ray flux in the local interstellar medium, after propaga- 
tion in the Galaxy from multiple sources. In the present work, 
the amount of particles injected into diffusive shock accelera- 
tion is taken as a constant fraction of the material crossing the 
shock. We used different injection fractions for protons and elec- 
trons, T]inj, p and rii„j e , and these were treated as free parameters. 
Particles enter the acceleration process with a fixed momentum, 
which we chose to be p mi =1 Me V/c. The exact value has no 
real importance because of a degeneracy in our model such that 
the level of non-thermal energy is proportional to the product 
TjinjPmj. We emphasise that constant injection parameters over 
space and time is a simplification and that they are very likely 
depen dent on the shock conditions dBlasi et al.ll2005t IVolk et all 
l2003h . 

The non-thermal populations produced by DSA simultane- 
ously undergo momentum losses due to adiabatic expansion 
of the zones where particles are confined, and radiative losses 
through a variety of processes. Electrons are subject to syn- 
chrotron and relativistic Bremsstrahlung emission in the dense 
and magnetised shocked material, and inverse-Compton scat- 
tering in the intense radiation field of the red giant and nova; 
protons experience inelastic collisions with gas particles in the 
downstream medium, thereby creating unstable mesons that de- 
cay in y-rays (see Sect. 13. 3I >. The physical conditions used to 
compute radiative losses at each time are those of the immediate 
post-shock medium. For zone A, this can be justified by the fact 
that particles in the acceleration process spend more time down- 
stream than upstream; for zone B, this is a simplification since 
successive generations of advected particles experience in prin- 
ciple a continuum of conditions as the shocked gas flows away 
from the shock front. 

The evolution of the non-thermal particle distribution N(p) 
in each zone is computed through the following general equation 



(with subscripts A or B for zone A or B): 

8N a ,b 



dt 



= T" KPa.bNa.b) + Qa.b 

op t a ,b 



(9) 



where p is the momentum gain/loss rate, t the characteristic es- 
cape time, and Q the source term. In the acceleration zone, these 
quantities have the following expressions: 



A \ dt ) D SA \ dt ladiab \ dt had 



T A = 



' DSA 
q(\ +q)D{p) 



riinj dM m 

Qa = -T—Oip-Pinj) 

pjfip at 



(10) 

(11) 
(12) 



where m p and p. are the proton mass and the mean molecular 
mass, respectively. Particle momentum evolves under the com- 
peting effects of diffusive shock acceleration, adiabatic expan- 
sion, and radiative losses (the subscript rad refers to one or a 
combination of the processes introduced above, depending on 
whether the particle is a proton or an electron). The escape time 
is computed as the ratio of the average cycle time by the escape 
probability. The source term is a Dirac at the momentum pi„j. In 
the cooling zone, the governing equation is: 



T B = OO 

a - Na 

<Jb - — 

TA 



(13) 
(14) 
(15) 



The source term is now the escape from zone A and there is no 
acceleration there. Once in zone B, particles can only accumu- 
late and lose momentum. 

The environmental conditions just ahead of the shock wave 
are computed at each time from the following assumptions: a 
spherically-symmetric wind density profile around the RG and 
possibly a circumstellar density enhancement around the WD, as 
presented in Sect. 13. II with a composition of 90% hydrogen and 
10% helium in particle number; a magnetic field assumed to be 
in equipartition with the thermal energy density of the upstream 
material, taken to be highly turbulent and to be at a constant 
temperature T w ; the spherically-symmetric blackbody radiation 
fields of the red giant and of the burning envelope of the WD. 

3.3. Radiation 

The model described above was used to compute the emissions 
associated with the shock. The main objective of this work is 
to account for the non-thermal y-ray emission detected in the 
Fenra/LAT range, which can result from inverse-Compton scat- 
tering, Bremsstrahlung, and hadronic interactions. We also com- 
puted the synchrotron emission because radio observations in 
the ~ 1-10 GHz band could be an additional constraint on non- 
thermal processes. Last, we evaluated the thermal radiation from 
the shocked material, as X-ray observations can constrain some 
parameters of the binary system. 

Radiation fields: The red giant star light was modelled as a 
Planck spectral distribution with a temperature of 7#g=2600K 
and a radius Rrg = 500 R Q , corresponding to a total luminosity 
Lrg = 1O 4 L . With a peak V-band emission nearly 5 magni- 
tudes brighter than the usual level of emission from the system 
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(IMunari et al 1201 lh . li ght from the nova eruption can provide a 
significant contribution to the ambient photon field. An envelope 
forms rapidly around the white dwarf after nuclear burning is ini- 
tiated. Peak visual light corresponds roughly to m aximum exten - 
sion of the photosphere, typically 10 12 - 10 13 cm (Warner 2003). 
Thereafter, the envelope shrinks while maintaining a nearly con- 
stant bolometric luminosity L„ ova close to the Eddington lumi- 
nosity before mass loss and the end of nuclear burning terminate 
the process. Detailed models of fast nova eruptions on massive 
white dwarfs dKato & Hachisul 19941) show that peak visual light 
is reached quickly after ignition, at which point the radius of 
the photosphere around the white dwarf is R„ ova — 5 x 10 12 cm 
and has a temperature T lmva = 10 4 K for a bolometric lumi- 
nosity L nova = 2 x 10 38 ergs -1 = 5 x 10 4 L Q . We found that 
a very good fit to the initial 50 days of the observed V band 
lightcurve is obtained by assuming blackbody emission and a 
photosphere shrinking as R nova — 5 x 10 12 (?/0.1 day)~'/ 2 cm 
(keeping L nova constant, hence T lmva oc f 1 ^ 4 ), for an adopted dis- 
tance of 2.7 kpc. This simple phenomenological approach allows 
to calculate the nova radiation field at all points and times. The 
red giant and nova luminosities are only a factor 5 different. The 
nova field therefore dominates during the early evolution as it is 
closest to the shock. Radiation from the shock-heated material 
was neglected because its observed luminosity (in X-rays) is a 
few 10 34 ergs~', so that the radiation density is no higher than 
10" 4 that of the nova for R s > 5 x 10 cm; in addition, its con- 
tribution to the scattered spectrum would be suppressed by the 
Klein-Nishina cutoff of the inverse Compton cross-section. 

Non-thermal emission: For the leptonic radiation processes, 
synchrotron and inverse-Compton radiation of high-energy elec- 
trons were computed usi ng the standard formula for energy loss 
rate and spectra from (Blu menthal & Gouldlll970T) . As shown 
below, relativistic Bremsstrahlung is negligible and its contribu- 
tion to the spectral energy distribution was not computed. The 
synchrotron emission spectrum was computed with the approxi- 
mation of a 90° pitch angle (which overestimates the total emis- 
sion power by 3/2). For inverse-Compton scattering on the red 
giant light, the transition to the Klein-Nishina regime occurs at 
a particle energy ~400 GeV. For the nova light, the transition is 
at ~ 100 GeV at early times, when nova emission is most rel- 
evant. As shown below, accelerated electrons will reach ener- 
gies of a few 10 GeV at most, so inverse-Compton emission in 
the Fermi/LAT range is thus in the Thomson regime. Inverse 
Compton emission and losses were calculated in the isotropic 
approximation for most runs, but we also investigated the in- 
fluence of anisotropy of the seed photon field on the spectrum 
se en by differen t obser vers (see Sect. 15.41) . following the method 
of iDubus et alJ (2008). In this case, the anisotropic inverse- 
Compton calculation includes the finite size of the red giant 
or nova photosphere and requires a discretization in azimuth 
to take into account the change in the scattering geometry. For 
the hadronic radiation processes, pion decay emission was com- 
puted following the parameterization of lKamae et alJ d2006l) . We 
did not use any nuclear enhancement factor to account for the 
contribution of nuclei heavier than protons (which would in- 
crease the emission le vel by a factor of about 2 for a solar metal- 
licity; see Morill2009l) . 

Thermal emission: Thermal Bremsstrahlung of the mate- 
rial heated by the shock is responsible for the X-ray emis- 
sion. The typical mean density of the swept-up material down- 
stream of the shock is computed assuming the shell width 
AR S is 10% of the shock radius. This v alue appears reason - 
able based on the numerical simulations o f lSchureetalJ feOlO). 
The thermal Bremsstrahlung spectrum was calculated follow- 



ing lRvbicki & Lightman (Il986l) . using as post-shock density the 
mean value defined above and as post-shock temperature 

T s = ^rum p V* * 10 8 ( — -\ K (16) 

or kT s » 8. 6 keV. The whole a pproach is very similar to what 
was done in iNelson et al.l(l2012l) . Yet, it suffers from various lim- 
itations such as the fixed shell width, or the absence of heat con- 
duction and gas mixing. Below, the comparison with compre- 
hensive hydrodynamical simulations of V407 Cyg will reveal 
these shortcomings and tell to what extent our thermal model 
can be relied on. 

Absorption: We have taken into account synchrotron self- 
absorption and thermal free-free absorption within the shell 
(again assuming AR S /R S = 0.1). Absorption within the red giant 
wind is discussed in Sect. 16.21 For simplicity, we have reduced 
the radiative transfer in the shell to a plane-parallel approxima- 
tion. Internal absorption is typically small except at low frequen- 
cies < 1 GHz and therefore affects only the synchrotron emis- 
sion. The flux after absorption for each position along a given 
propagation angle is computed as 

^2 ^^sync^Tfree^ 

/abs = (/sync + /free) ~ (17) 

fsync Tfree 

where the opacities were compu ted from the local co n dition s 
in the shock element following Rvbicki & Lightman (1986). 
Absorption of y-rays due to pair production was neglected be- 
cause the density of X-ray photons from the shock (the only 
photons that can interact with GeV photons to create pairs) is 
too low. Assuming that the shocked material emission is all in 
1 keV photons, then their density is nx < 10 7 phcrrT 3 , leading 
to a pair production opacity r ~ crjnxRs ^ 10~ 6 . We have also 
ignored X-ray absorption of the shock emission. The typical col- 
umn density for the stellar wind integrated down to the surface 
of the red giant is Nh ~ 10 23 cmT 2 , which leads to absorption 
below 3 keV. Thermal Bremsstrahlung emission is harder at the 
beginning and the column density sampled becomes smaller as 
the shock progresses outwards such that absorption of the shock 
emission after 20 days above a keV or so, in the range where 
we compare our results to the Swift/XRT observations, can be 
neglected. 

3.4. Time scales 

We give below some estimates for the typical time scales of the 
different processes involved to facilitate the interpretation of the 
results presented in the following sections. We used as average 
environmental conditions those found at the WD position in the 
base case scenario defined in the next section. The particle den- 
sity, radiation energy densities, and equipartition magnetic field 
strength upstream of the shock are respectively 

n gas = 3 x 10 7 cnT 3 (18) 

u RC = 4 x 10~ 3 erg cm" 3 (19) 

u nova = 20 erg crrT 3 (20) 

B. S = 10- 2 G (21) 

where the nova radiation energy density is calculated at the peak 
assuming the shock is located at the photosphere. The nova ra- 
diation field at the shock location decreases rapidly as it moves 
away (initially u„ ova oc r 2 ), so that the mean u nova and mrg typ- 
ically become comparable within a week. The particle density 
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used in energy loss calculations is the downstream one and is 
thus a factor of 4 larger than the above value, at about 10 8 cm -3 , 
because of shock compression. 

To begin with, the timescale to accelerate a particle to an 
energy £=10 GeV is 



fDSA *°^(lo!ev) 



0.01 G\ /3000 km s" 



B s 



V s 



day (22) 



for the typical initial shock velocity and magnetic field intensity 
experienced by particles over the first days or weeks (depending 
on propagation angle). At the lowest energies, Coulomb losses 
could compete with DSA owing to the relatively high densities. 
Coulomb losses actually dominate over ot her loss proces ses be- 
low ~100MeV and occur on a time scale (lLongairll994 

Ae /10 8 cm- 3 \/ E \ J 

fc — * 45 (^H(lb^) day (23) 

This is much longer than the acceleration timescale, so we have 
neglected its influence. Similarly, adiabatic decompression due 
to shell expansion could compete with DSA especially in the 
early stages. The corresponding time scale for a typical shock 
radius and velocity is 



^adiab ~ 3 



V5AU/1 



3000 km s" 

V s 



day 



(24) 



which is more than one order of magnitude above the time scale 
for acceleration (and even more for lower energies). Particle ac- 
celeration can therefore proceed up to higher energies, where it 
competes with other energy loss processes. Focusing on elec- 
trons first, the synchrotron loss timescale is 



5000 



0.01 G \ / 10GeV \ 
Bs ) \ E } 



day 



(25) 



In the case of significant magnetic field amplification, the value 
for B s would have to be augmented but we will see below that 
strong amplification of the mean magnetic field is not required. 
The energy density in stellar radiation is about 1000 times the 
magnetic energy density, so inverse-Compton can be expected 
to dominate synchrotron losses and be characterised by a far 
smaller time scale 



fic 



4 x 10 3 erg cm 



-3 



10GeV\ 
—E~ day - 



(26) 



where inverse-Compton losses were computed in the Thomson 
approximation. The last process responsible for electron energy 
loss is relativistic Bremsstrahlung, which occurs on a time scale 



dLongairil 19941) 



/10 8 cnr 3 \ 
f Br ~ 300 | | day 



(27) 



which is far longer than that of inverse-Compton. Now look- 
ing at protons, energy losses through inelastic collisions with the 
shocked gas proceed over a characteristic time 



f pp » 260 



10 8 cm" 



l gas 



day 



(28) 



where we have used a threshold for the reaction at 280 MeV, a 
constant inelasticity parameter k pp = 0.5, and a constant cross- 
section cTp P = 30 mbarn (which is a reasonable approximation 



for proton energies below 1 TeV, see lKelner et alJ [2006). 

Overall, for the adopted average conditions, we can expect 
particle acceleration to proceed up to ultra-relativistic energies. 
For electrons, the initial acceleration is rapidly stuck at 1 GeV 
by inverse-Compton losses on the intense nova radiation field 
at the early shock position; it later continues up to energies of 
order 100 GeV, when electron acceleration is limited by inverse- 
Compton losses on the lower-density radiation fields found at 
larger distances from the binary. For protons, acceleration could 
continue up to energies of order 10 TeV, but we will see be- 
low that it stops before as a result of shock slowing-down. 
Electrons will therefore very likely be loss-limited and protons 
age-limited. In terms of radiation in the 100 MeV- 1 GeV range 
(which results from 10 GeV particles, either through inverse- 
Compton or hadronic interactions), non-thermal electrons are at 
least ~100 times more efficient in radiating away their energy 
than non-thermal protons. Yet, if the latter are at least 100 times 
more numerous than the former, they will dominate the y-ray 
emission. 

For the thermal emission, the cooling timescale of the shock- 
heated wind is 



tfs ~ 170 



(io 8 k) 



1/2 / 10 



day 



(29) 



using a Gaunt factor equal to 1, and T s as defined previously. 
The cooling timescale of the shocked material is long, which 
justifies the assumption of Sedov-Taylor expansion in the first 
couple of months, provided most of the shock kinetic energy 
is indeed channele d into heating up the gas thermally (see also 
lNelsonetalJl2012h . 



3.5. Computation 

The evolution of the whole system was computed for a grid of 
propagation angles over successive time steps, using an angular 
step of A9 = 7r/100rad and a time step of Af = 0.1 d. The non- 
thermal particle populations were followed on a momentum grid 
running from /?„„■„ = 10 -1 to p max = 10 8 with a logarithmic step 
of 0.01 (in units of MeV/c). 

For a given set of parameters, the evolution is computed 
along each direction 9, starting at the position of the WD sur- 
face with the initial shock velocity. At each time step for a given 
propagation angle: the shock properties (radius, velocity, swept- 
up mass, dynamical stage) for the current shock element are up- 
dated; the environmental conditions are determined (gas and ra- 
diation density, magnetic field strength); the non-thermal parti- 
cle distributions are computed according to Eq.|9] first in zone A 
and then in zone B. All along, we check that the current shock 
element does not crash into the RG by a simple criterion on posi- 
tion. If that is the case, the calculation is stopped for the current 
direction and the associated particle population is removed from 
the total. We did not handle the interaction of the shock with the 
RG surface or the convergence of the shock on the rear side of 
the star, so the shadow cone behind the RG is a dead zone in our 
model. Once the hydrodynamical and non-thermal evolution is 
computed, the corresponding radiative outputs are determined in 
a post-processing stage. 

For a given layout of the V407 Cyg binary system, there 
are three free parameters: the diffusion efficiency, and the injec- 
tion fractions for protons and electrons. Since we neglected the 
non-linear aspects of the acceleration process, the calculation is 
linear and the injection fractions are just normalisation factors, 
tuned to match the y-ray observations. In any case, they are con- 
strained by energetic considerations: the energy channelled into 
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10° 10 2 10" 10 6 10 s 10° 1 2 10* 10 6 10 s 

Momentum (MeV/c) Momentum (MeV/c) 

Fig. 2. Time evolution of the proton and electron distributions (top and bottom row, respectively) in the acceleration zone and cooling 
zone (left and right column, respectively). The results correspond to the base case scenario with energy losses, a propagation angle 
of 47.3° from the WD-RG axis, and injection fractions of r\mj,p = 10~ 3 and r\i n j, e = 1CT 5 . Curves from black to red show successive 
time steps of 0.5 day. The dashed curve in each plot marks the ED to ST transition. 



non-thermal components should not exceed a fraction of the to- 
tal nova kinetic energy. As we will see, the diffusion efficiency 
provides control on the maximum particle energy, and thus on 
the y-ray spectral maximum. The value eventually required for 
this parameter provides information about the particle scattering 
regime and possible magnetic field amplification. 

We present below the outcome of three different calculations, 
in order to emphasise how thermal and non-thermal emissions 
depend on the system properties. The first calculation is based 
on typical para meters for the V407 Cyg system, quite similar to 
those adopted in lAbdo et al.l (2010), with a surrounding medium 
consisting only of the RG wind. This base case scenario will be 
useful in illustrating the main features of particle acceleration 
in such a system but cannot account for several observations. A 
second calculation is then presented, where the binary system 
and acceleration parameters were optimised in order to reduce 
the discrepancy with the experimental data, but still keeping a 
wind density profile for the ambient medium. This scenario also 
has major shortcomings and thus calls for improvement of the 
model. In a last calculation, a density enhancement was added 
to the simple wind profile and yielded a far better fit to most ob- 
servations. These calculations will be referred to as Run 1, Run 
2, Run 3, respectively, and their corresponding sets of param- 
eters are summarized in Table [T] In terms of particle accelera- 
tion, we aimed at accounting for the ferm//LAT data with the 
smallest amount of non-thermal energy, in order to remain close 



to the limit of test particle approximation, and with the smallest 
electron-to-proton ratio at injection K ep = 7/,- n / ;g / 7/,- n ; - >p , in order to 
be in agreement with the values quoted in Sect. 13.21 In the base 
case scenario, for illustration, particle acceleration was assumed 
to occur in the Bohm limit in the non-amplified equipartition 
magnetic field (the field strength is already quite high, typically 
about 10 4 times the local interstellar value), and the electron-to- 
proton ratio was set to 1%. 

4. Non-thermal particle populations 

In this section, we focus on the non-thermal particle populations 
that can be produced by V407 Cyg in the base case scenario 
(as defined in Sect. 13.5b . We present their evolution in time and 
space over the first 40 days after the nova outburst. We first dis- 
cuss the properties of the particle distributions especially their 
maximum energy and the impact of energy losses, then the ef- 
fects of the asphericity of the phenomenon on particle accelera- 
tion, and last the global energetics. 

4. 1 . Distribution properties 

In Fig. [2] the time evolution the non-thermal particle popula- 
tions in the acceleration and cooling zones is plotted, for a prop- 
agation angle of about 45°. The distributions obtained for other 
directions are very similar, except the ones for low angles (prop- 
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Table 1. Summary of the different calculations presented in the 
paper. 



Run 


Binary system parameters 


Acceleration parameters 




M ej = 1(T 6 M 


77,,,^ = 3 x 10- 3 


i 
1 


v e j— 30U0 Km s 


?7mj, e = 3 X 10 




d orb = 10AU 


# cp = 0.01 




M RG = 5 x 1(T 7 M yH 






Rwd= 0.01 R 






/J fiG =500R Q 






L m = 10000 L Q 






y - 70 km « _I 






T - 700 K 

± yrf / \J\J 1^ 






PCDE — u v i_-T, ; 






M e ;= 2 x 10~ 6 M Q 


r}inj, P = 6 x 10~ 3 


2 


M SG = 10- 7 M o yr-' 


r, inM = 6 x 10- 4 




rf ori ,= 6AU 


K ep = 0,l 




p CO£ = (no CDE) 






M ej = 2 x 10~ 6 M 


77,,,^ = 5 x 10~ 3 


3 


M RG = 5 x 10- 8 M o yr-' 


77,,,^ = 3 x 10- 4 




rf OT fc= 10AU 


= 0.06 




Pcde = 10 s cm" 3 






Icde = 10 AU, W =5AU 





Note to the table: CDE stands for Circumstellar Density Enhancement 
(see Sect. l3.lt . For Run 2 and Run 3, the parameters not listed have the 
same values as in Run 1 . 



agation towards the RG), whose evolution is prematurely halted 
by the collision with the star. 

We first discuss the case of accelerated protons, which are 
almost only exposed to acceleration and downstream advec- 
tion, hadronic interactions and adiabatic losses being not effi- 
cient enough. Therefore, they approximately illustrate the case 
of acceleration without energy losses. In the acceleration zone, 
relativistic particles soon follow a power-law distribution in mo- 
mentum with a slope s = — 1 and a cutoff at an increasing mo- 
mentum. In the cooling zone, relativistic particles also follow 
a power-law distribution with cutoff, but steeper with s = -2 
because of the energy dependence of the advection time, which 
favours the escape of lower-energy particles (see Eq. [TTlPl The 
cutoff in both distributions is set by the finite acceleration time, 
hence its progressive shift to higher and higher momenta. This 
is especially true for acceleration during the ED-stage, where 
the acceleration rate remains almost constant because of a con- 
stant shock velocity. The progression to higher energies seems 
to stall because the acceleration rate is independent on the par- 
ticle momentum (see 13.2b . and so a given logarithmic increase 
in energy takes longer to be achieved at higher momenta. In ad- 
dition, for those parts of the shock front that enter in ST-stage, 
the acceleration rate decreases because of shock slowing-down, 
and the cutoff momentum hardly evolves anymore. It is however 
during the ST-stage that most particles are accelerated. Overall, 
non-thermal protons are almost age-limited. In the base case sce- 
nario, the maximum particle energy is slightly short of 100 GeV 



3 The proton distributions for zone A exhibit a break at about 1 GeV/c 
that marks the transition to the relativistic regime. This feature arises 
because of the /3 factor in the diffusion coefficient D(p), which implies 
a very efficient acceleration of mildly-relativistic particles. In the same 
time, this means an easier downstream advection of the same particles, 
so both effects compensate and the break does not appear in the proton 
distributions for zone B (seeEas.l8land[TTI 



and non-thermal emission in V407 Cyg 
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Fig. 3. Spatial maps of the maximum energy attained by protons 
(top) and electrons (bottom) in the base case scenario. The color 
scaling is linear between the highest value (in black) and values 
that are less than 10 times smaller (in white). Note that the max- 
ima are different for protons and electrons (see text). The WD 
is in (x,y)-(0,0) and the RG in (x,y)=(10,0). Dotted white lines 
mark successive positions of the shock front with a 2-day time 
step, and the solid white line marks the ED-ST transition. 



after just 1 day, and it increases by about one order of magnitude 
to almost reach 1 TeV when the shock starts to slow down. After 
that, the maximum particle energy remains almost constant up 
to at least day 40. 

The situation is somewhat different for accelerated electrons. 
Their distribution is strongly modified by inverse-Compton 
losses, as illustrated by Fig. [2] In the acceleration zone, elec- 
trons are initially stuck at GeV energies because of the very 
strong inverse-Compton losses on the nova light (if they were 
exposed to the RG light only, they would move up to about 
50 GeV in just one day). The competition of acceleration and 
inverse-Compton results in a pile-up at the high end of the distri- 
bution. As the shock moves away from the nova photosphere, the 
inverse-Compton loss rate drops and the particles can be pushed 
to higher energies, eventually reaching 30 GeV at the end of the 
ED-stage. In the ST-stage, the acceleration rate decreases and 
inverse-Compton losses prevail, which causes a progressive shift 
of the pile-up to lower energies. This behaviour is reflected by 
the distribution in the cooling zone, which in addition is altered 
at the high end by inverse-Compton losses and exhibits a char- 
acteristic steepening from s = -2 to s = -3. 

4.2. Effect of asphehcity 

The aspherical shock wave and anisotropic binary environment 
have an impact on accelerated proton and electron energies and 
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Fig. 4. Spatial map of the non-thermal energy stored in acceler- 
ated protons in the base case scenario (the plot for electrons is 
quite similar). Color coding and graph elements are described in 

Kg. m 



on how the non-thermal energy is distributed over the shock 
front. This is illustrated in Figs.[3]and[4] 

For protons, there is a mild dependence of the maximum par- 
ticle energy on the direction of shock propagation. The transition 
from ED to ST-stage occurs later with increasing propagation 
angle, but this has a limited impact because of the equiparti- 
tion assumption for the magnetic field: a lower average density 
along a given propagation direction ensures a longer ED-stage, 
but implies a smaller magnetic field; eventually, a lower accel- 
eration rate operates for a longer time, giving about the same 
maximum energy. As a result, the maximum proton energies are 
quite uniform after a few days and remain almost unchanged 
subsequently. The highest particle energies are reached for prop- 
agation towards the densest regions because of stronger equipar- 
tition magnetic fields. 

For electrons, which are strongly affected by energy losses, 
the picture is different. For shock elements propagating towards 
a denser radiation field, inverse-Compton losses limit the maxi- 
mum particle energy; the highest energies are therefore reached 
for propagation away from the RG. As a result, the maximum 
electron energy is achieved over a different domain of the shock 
trajectory and it is less uniform at a given time than for protons. 
Another difference is that the maximum energies will decrease 
past the ED-ST transition whatever the propagation direction, in 
contrast to protons. 

Because of the relatively limited variation in maximum par- 
ticle energy over the shock surface, especially for protons, the 
non-thermal energy distribution over the shock front is essen- 
tially set by the number of accelerated particles, hence the 
amount of swept-us mass. For a blast wave propagating in a red 
giant wind only, shock elements moving along angles between 
about 10° and 90° are those that dominate the production of non- 
thermal particles (see Fig.|4|. 

4.3. Global energetics 

Particle acceleration channels a fraction of the nova kinetic en- 
ergy to non-thermal components. This energy first goes into 
non-thermal particles, which transfer it partially to radiation and 
magnetic fields. In the following, we will define at each time the 
fraction of the nova energy that is in the form of non-thermal 
particles as acceleration efficiency, and refer to the fraction of 
non-thermal particle energy that has gone into non-thermal radi- 



ation as radiation efficiency. Then, the non-thermal efficiency is 
defined as the fraction of the nova kinetic energy channelled into 
non-thermal components, particles or radiation. In the energy 
budget, we neglect the magnetic energy drained from accelerated 
particles by wave amplification processes (this can extract <10% 
of the incoming flow energy for parallel, non-relativistic s hocks 
with Mach number M A > 30: lGargate & Spitkovskvl2012l) . 

As emphasised previously, the non-linear aspects of the par- 
ticle acceleration process were not considered in our model. The 
main consequence is that the energy and flux levels of non- 
thermal components are proportional to the injection fractions 
(for a given injection momentum). For the values iji„j, p = 3xlO~ 3 
and rji„j e — 3 x 1CT 5 required to match the observed y-ray flux 
(see below), the acceleration efficiencies at day 40 for protons 
and electrons are about 100% and 2%, respectiveljQ. Such injec- 
tion fractions in the base case scenario therefore seem to be ex- 
cluded by energetic considerations. As we will see below, how- 
ever, the correspondence between injection fraction and acceler- 
ation efficiency depends on the environmental conditions. 

Last, it is interesting to compare the content of both zones. 
The cooling zone contains more particles and energy than the 
acceleration zone, but the situation is different for electrons and 
protons. After a few days, there is about 4 times (respectively 30 
times) more non-thermal energy in the cooling zone than in the 
acceleration zone for protons (respectively electrons). This dif- 
ference is not uniform across the energy range due to the differ- 
ent spectral distributions, as can be seen from Fig. [2] The cooling 
zone contains about a few 100 times more ~1 GeV particles than 
the acceleration zone, but this contrast diminishes with energy 
and the contents become comparable at ~ 100 GeV. This is so 
because higher-energy particles take longer to be produced and 
have a longer escape time. And since protons can reach energies 
that are 10-20 times higher than that of electrons, they can store 
more energy in the acceleration zone than electrons. 

5. Gamma-ray emission 

In this section, we discuss the y-ray emission associated with the 
V407 Cyg outburst. We first show the results obtained in the base 
case scenario, where the shock propagates in the RG wind with 
typical parameters for the binary system and acceleration pro- 
cess. This illustrates several properties of the high-energy emis- 
sion from such a system but, as we will see, it cannot account 
for the observations of V407 Cyg. We then present an attempt to 
optimize the binary system and acceleration parameters in order 
to reduce the discrepancy with the experimental data, and from 
that conclude that the observations cannot be reproduced for a 
shock propagating in a wind only. In a third part, we improve 
our model by considering a more realistic ambient medium that 
includes a density enhancement on top of the wind profile. 

5.1. Shock propagating in a wind: base case scenario 

Based on the non-thermal particle distributions presented in 
Sect. H] we computed the y-ray emission of the V407 Cyg out- 
burst due to inverse-Compton scattering of the electrons and 
hadronic interactions of the protons. The results are presented in 
Fig. [6] for both populations simultaneously. The injection frac- 

4 The acceleration efficiency for electrons is not 100 times less 
than the acceleration efficiency for protons, as could have been ex- 
pected from the ratio of the injection fractions, because at low energies 
< 10 GeV electrons gain more kinetic energy than protons for the same 
momentum gain. 
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Fig. 5. Spatial maps of the >100MeV y-ray emission frompion 
decay (top) and inverse-Compton (bottom) in the base case sce- 
nario. Color coding and graph elements are described in Fig. [3] 



tions were set at r\mj tP — 3 x 10~ 3 and r\ m j, e — 3 x 10~ 5 , in 
order to match the level of the y-ray emission. We discuss be- 
low the properties of these emissions and compare them to the 
FermifLAI data. We emphasise that this discussion is really 
specific to the injection fractions and binary system parameters 
adopted in the base case scenario. 

Spectrum: Below 150MeV, the emission is contributed in 
similar proportions by inverse-Compton on the nova light and 
RG light. Above 150MeV, the emission associated with neutral 
pion decay dominates over that from inverse-Compton. The to- 
tal time-averaged spectrum exhibits a broad maximum at about 
a few GeV and exceeds the highest-energy points from the 
Fermi/LAT measurement. This results from a non-thermal par- 
ticle population extending too high in energy. A lower momen- 
tum cutoff in the particle distribution would sharpen and shift 
the spectral profile to lower energies, making it more consistent 
with the observations. This could be obtained with an enhanced 
diffusion efficiency, as shown below. 

Light curve: The time-averaged photon flux is at about the 
expected level, but the predicted light curve clearly differs from 
the observed one. The early rise is contributed almost exclu- 
sively by inverse-Compton on the nova light, but at too slow 
a rate. The flux maximum is reached at about day 6-7, while 
y-ray observations indicate a maximum at day 3-4. Moreover, 
the light curve does not fall as rapidly as observed: at day 40, 
the predicted photon flux has dropped by a factor of ^5 from 
its peak value, whereas y-ray observations showed a decrease 
by more than an order of magnitude (the collective upper limit 
on the emission from day 19 to 33 is 0.8 xl0~ 7 phctrT 2 s _1 , not 
shown to not overload the plots). The emission peak occurs as 
the shock moves through the regions of highest matter and radi- 
ation density, which in the base case scenario are found close to 



the RG surface. The subsequent shallow decrease occurs because 
as the shock moves away from the binary, the steady increase of 
the number of particles with energies of ~ 1 0- 1 00 GeV (those re- 
sponsible for y-ray emission in the ~l-10GeV range) partially 
compensates the decreasing matter and radiation density. 

Locus: Figures [3] and |4] show that the maximum momentum 
and energy budget of accelerated particles vary modestly over a 
sizable fraction of the shock front. In contrast, the high-energy 
emission is contributed by a limited part of the blast wave. This is 
illustrated by Fig.[3]for the inverse-Compton and pion decay ra- 
diation processes. Both emissions are concentrated in the vicin- 
ity of the star, where the number of accelerated particles and 
the density of target gas and photons are the largest. Due to its 
higher radiative efficiency, however, the leptonic emission arises 
from a more extended domain. This explains why the pion decay 
light curve drops slightly faster than the inverse-Compton one. 
This shows once again that in the base case, the flux maximum 
is reached approximately when the shock reaches the RG. 

Energetics: With the injection fractions adopted, the total ac- 
celeration efficiency at day 40 in the base case is about 100 %, 
with protons largely dominating the non-thermal energy budget. 
The non-thermal energy output to radiation is about the same for 
protons and electrons, despite the former being far more numer- 
ous than the latter. This is due to the higher radiation efficiency of 
high-energy electrons. The bottom plot of Fig.|6]shows that pro- 
tons radiate inefficiently even in the dense RG wind and retain 
most of their non-thermal energy (only =* 1 % has been converted 
in radiation after 40 days), while electrons lose their non-thermal 
energy quite rapidly (about ^20% have been converted in radia- 
tion after 40 days). 

Overall, with the base case hypotheses, our model gives a y- 
ray emission that is at approximately the observed level, and the 
y-ray energies are about the measured ones for Bohm diffusion 
in the wind equipartition magnetic field. The predicted emission, 
however, reaches its flux maximum too late and does not decay 
rapidly enough. More problematic are the non-thermal energet- 
ics: at the time of the predicted emission maximum the accelera- 
tion efficiency is a reasonable ^20%, but the continuation of the 
acceleration process leads at day 40 to a total non-thermal en- 
ergy in excess of the initial nova kinetic energy. Potential solu- 
tions to this energy problem could be a larger ejecta mass, which 
means a larger kinetic energy reservoir, or a higher electron-to- 
proton ratio at injection K ep , where the proton injection would 
be reduced to not put a strain on the energy budget, while the 
electron injection would be increased to provide enough y-rays. 
Yet, a more prohibitive shortcoming is the associated thermal 
emission in X-rays (not shown), which is about 30 times above 
the luminosity maximum inferred from Swift/XKY observations. 
This means that the medium density is too high and that too 
much material is swept-up by the blast wave. If the environment 
just consists of the RG wind, the stellar mass-loss rate should be 
Mrc ^ 10~ 7 M yr~' . In the next subsection, we investigate how 
changing the main parameters of the system impact the emission 
properties, and we search for a combination of plausible values 
that could lead to a better agreement with the observations. 

5.2. Shock propagating in a wind: parameter dependence 

The initial shock speed can be considered to be established with 
relatively high accuracy from emission line observations, but 
most other parameters have uncertainties reaching about an or- 
der of magnitude in some cases. Below, we examine their ef- 
fects on particle acceleration and non-thermal emission, chang- 
ing only one parameter at a time from the base case scenario. 
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Orbital separation: Decreasing the orbital separation means 
a shock starting its propagation in higher gas densities. The 
ED-stage is shorter and the shock velocity drops faster in ST- 
stage, but in terms of particle acceleration this is compensated by 
stronger magnetic fields, hence smaller diffusion coefficient and 
higher acceleration rate. Increasing the orbital separation leads 
to exactly the opposite. We tested three values R or t,-5, 10,15 AU, 
and in all cases about the same maximum energy is reached by 
non-thermal particles, and about the same number of particles 
have been accelerated at day 40. As to the y-ray emission, the 
main difference is in the light curve: for the smaller orbital sepa- 
ration case, the maximum particle energy is reached sooner and 
closer from the star; as a consequence, the light curve reaches 
its maximum sooner and the peak is more pronounced. A larger 
orbital separation leads to the opposite. Apart from that, all three 
values of R or \, yield about the same average spectrum, and a rel- 
atively flat light curve past maximum. 

Mass-loss rate: The density encountered by the shock de- 
pends linearly on the mass-loss rate. We tested several values 
around lO _7 M yr _1 . A lower mass-loss rate means a longer 
ED-stage but a lower magnetic field, and this eventually gives 
a similar maximum particle energy. Yet, with a lower mass-loss 
rate, higher-energy particles are produced at a greater distance 
from the binary system, and are taken further out at a greater ve- 
locity. This means that they are exposed to lower radiation and 
gas densities and for a shorter time, which implies a lower ra- 
diation efficiency. The low density affects in particular protons, 
whose relative contribution to the y-ray emission can fall below 
that of electrons. In addition, a lower wind density for a same 
injection fraction also results in fewer particles being acceler- 
ated and a lower acceleration efficiency at a given time. In such 
conditions, reproducing a given non-thermal emission level at 
early times requires an increase of the injection fractions, and 
may lead to unrealistically high acceleration efficiencies at late 
times. In terms of y-ray emission, a lower mass-loss rate leads 
to a sharper and earlier light curve maximum, with a more rapid 
decline but a smaller flux overall (for a same injection fraction). 
A higher mass-loss rate leads to the oppositq^. 

Ejecta mass: The ejecta mass linearly sets the kinetic energy 
reservoir from which non-thermal components are powered. Too 
small a mass requires too high an acceleration efficiency to ac- 
count for a given level of y-ray emission, and we therefore con- 
sidered only ejecta masses greater than 10~ 6 M Q . A higher ejecta 
mass implies higher shock velocities on average, because the ST- 
stage slowing-down occurs later and on a longer time scale. At 
a given time and for the same injection fractions, more particles 
have been accelerated and they reach higher energies because ac- 
celeration at full shock speed proceeded on a longer time. In ad- 
dition, accelerated particles are advected away more rapidly, to- 
wards smaller gas and radiation densities. Despite these effects, 
the global acceleration efficiency is lower. The y-ray emission in 
the larger ejecta mass case features a more energetic spectrum, 
a higher luminosity, and a more pronounced light curve maxi- 
mum. This light curve maximum is however reached at about 
the same time, and the subsequent emission decrease remains 
nearly as flat. 

Diffusion coefficient: The base case calculation was based on 
the assumption of Bohm diffusion (in a non-amplified upstream 
field), which provides the maximum acceleration rate. This pre- 
scription can be relaxed through the parameter £ in Eq.[8] which 



s There is a certain degeneracy between the mass-loss rate and the 
wind velocity in the sense that the same density profile can be obtained 
through different combinations of both parameters. 



allows scaling up the scattering mean free path as a multiple 
of the Larmor radius. When energy losses can be neglected (as 
for protons), this results in a maximum particle energy that is 
lower by a factor but if energy losses play a major role (as 
for electrons), there is still a reduction but not by £ because the 
maximum energy is then set by a balance between acceleration 
and energy losses. In terms of energetics, increasing the diffu- 
sion coefficient decreases the total non-thermal particle energy. 
Proportionally, this mostly affects protons because they are age- 
limited while electrons are loss-limited. The consequences of a 
moderately higher diffusion coefficient on the y-ray emission are 
mainly a spectrum shifted to lower energies, while the lightcurve 
is left almost unchangecQ- 

5.3. Shock propagating in a wind: optimised scenario 

From the above discussion, the trend to follow in order to get a 
better agreement with observations would be (compared to the 
base case scenario): higher ejecta mass, smaller orbital sepa- 
ration, larger diffusion coefficient, and smaller mass-loss rate. 
These choices favour a rapid rise and decline of the light curve, 
a global acceleration efficiency down to moderate values, a lower 
y-ray spectral maximum, and thermal emission at the level ob- 
served in X-rays. We tested several possibilities along this trend 
and eventually identified an optimised scenario that we will refer 
to as Run 2, and whose parameters are listed in Table Q] 

The lower wind densities imposed by thermal emission con- 
straints have two direct implications: a reduced proton radiation 
efficiency, and fewer particles being accelerated for a given in- 
jection fraction. To obtain y-ray emission at a given level, these 
effects can be compensated by a higher proton injection fraction, 
but that means an enhanced consumption of the nova kinetic en- 
ergy. In Run 2, getting pion decay emission at the observed level 
requires r\mj,p — 3 x 10 _1 and a non-thermal efficiency above 
100% after just a few days. Overall, with such density condi- 
tions, the y-ray emission cannot be of hadronic origin. A leptonic 
origin can be considered, but only at the expense of a higher K ep : 
agreement with the data can be obtained for r\i n j, P — 6 x 10~ 3 , 
r linj,e — 6 x 10~ 4 , as illustrated in Fig. [7] 

The predicted time-averaged spectrum is fully consistent 
with the measured one, and this is obtained for particle scat- 
tering at the Bohm limit, with t, = 1 . The predicted light curve 
reaches its maximum at day 3-4, as observed, but falls short of 
the emission detected over the first two days. Actually, the rising 
part of the light curve is consistent with the data if one takes into 
account the 3-day uncertainty in the exact date of the outburst. 
The predicted emission, however, does not decline fast enough 
beyond the first two weeks and exceeds the collective upper limit 
at late times by a factor of 2. All this is achieved with a global 
acceleration efficiency of 30% at day 40, but still on the rise 
so that it reaches =s 50% at day 80. 

The above scenario yields about the best possible fit to the 
y-ray observations for a shock propagating in a wind density 
profile. The level of thermal X-ray emission matches that ob- 
served with Svw/f/XRT past day 20. The adopted ejected mass 
is in the rang e of acceptable va lues for a massive WD with 
Mwd ^ 1 .25 dYaron et al J 12005). The adopted orbital separa- 
ti on is smaller than the commonly-acce pted values: about 1 6 AU 
in lMunarietail (1 19901) and >12 AU in lshoreet afl (1201 lh . Yet, 



6 The effect of ( can also be mimicked through a modification of 
the wind temperature, since this partially defines the equipartition mag- 
netic field (the implications on the synchrotron emission are different 
though). 
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a smaller separation cannot be completely excluded considering 
the uncertainties on the orbital period (Ulisse Munari, private 
communication). Interestingly, some features of the predicted y- 
ray emission proved to be rather robust to order-of-magnitude 
changes of the model parameters. First, obtaining a spectrum in 
agreement with the observed one generally required a particle 
scattering close to the Bohm limit (1 < £ < 10). As will be 
discussed below, this is a useful element for the assessment of 
novae in symbiotic systems as a population of high-energy y-ray 
emitters. For V407 Cyg, another robust feature of all models that 
have been tested is that the predicted emission never exhibited 
the strong drop observed by Fermi/LAT about two weeks after 
the initial detection. This calls for the replacement or sophistica- 
tion of some of the assumptions adopted so far. In the following 
part, we show that taking into account a more realistic ambient 
medium improves the fit to the observations. 

5.4. Shock propagating in a wind with density enhancement 

As shown above, a simple spherically-symmetric wind profile 
prevents a quick rise of the y-ray emission, unless the orbital 
separation is small, and in any case results in a shallow decrease 
of the light curve past maximum. The latter shortcoming is con- 
nected with another problem, which is the steadily increasing 
non-thermal efficiency at late times. In this section, we mainly 
investigate the effect of alternative, more realistic density distri- 
butions. 

From hydrodynamic sim ulations of the symbiotic RS Oph 
system, IWalder et al.l d2008l) have shown that the quiescent ac- 
cretion phase between outbursts leads to a structured mass dis- 
tribution with enhancement in the orbital plane. This favours 
a poleward propagation of the shock when the system enters 
in outburst. This conclusion and its applicability t o the V407 
Cyg system were corroborated bv lOrlando & Drake! (1201 2l) . who 
found that a density enhancement on top of the wind density 
profile is needed to account for the X-ray light curve (see Sect. 

01. Due to the cylindrically-symmetric nature of our model, we 
could not explore the shock propagation and associated non- 
thermal emission in a realistic three-dimensional matter distribu- 
tion. The density enhancement was therefore modelled accord- 
ing to Eq. [2 which is a very good approximat ion at least for 
thermal X-ray emission cOrlando & Dra ke 2012). We placed the 
centre of this structure on the WD, in agreement with the simu- 
lations ofWalder et al.l (120081) that show a decreasing density in 
any direction from the WD (see their Fig. 3). 

The effect of such a circumstellar component can easily be 
anticipated. It will provide more material for acceleration and 
interaction over just a fraction of the shock front, in particular 
between the WD and the RG. This will allow for an earlier rise 
of the y-ray emission, and for a more pronounced drop when 
the shock exits the density enhancement. The mass accumulated 
around the WD is constrained by thermal emission observations, 
since this additional material swept-up by the shock will shine 
in X-rays. We tested several possibilities and found that a CDE 
with a typical size of 10 AU and a peak density of 10 8 cm -3 pro- 
vides a significantly better fit to the y-ray data. This new con- 
figuration will be referred to as Run 3, and the corresponding 
parameter set is given in Table Q] Compared to the previous Run 

2, the density enhancement allowed to come back to a larger or- 
bital separation and to further decrease the wind density through 
a reduced mass-loss rate. 

Even with a matter structure around the WD, however, repro- 
ducing the observed y-ray emission with hadronic interactions 
requires high proton injection fractions and still poses an ener- 



getic problem, with acceleration efficiencies being close to 100% 
for some parts of the shock (those propagating towards the RG). 
So once again only a leptonic option can be considered, still at 
the expense of a high K ep . Figure[8]shows the y-ray emission for 
Run 3 for j] in j 4 , = 5 x 10~ 3 , T] inJ:t , = 3 x 10~ 4 . 

The predicted y-ray spectrum required again particle scatter- 
ing close to the Bohm limit with f = 3 in order to match the 
measured one. The emission is mostly of leptonic origin and 
consists primarily of inverse-Compton scattering on the nova 
light. Hadronic emission dominates only above lOGeV but in 
the downturn part of the pion bump, so significant emission in 
the 100 GeV range and beyond is not expected. In terms of spec- 
tral variability over the first 15 days, the predicted spectrum ini- 
tially has a cutoff at ~100MeV and extends up to ~1 GeV in 3 
days. The spectral shape is then relatively stable over the follow- 
ing 10 days. At later times, 30-40 days, the spectrum is harder 
because the electron distribution flattens at high energy, as the 
inverse-Compton loss rate drops over time (see Fig.|2]i. Overall, 
our result ar e consistent with t he absence of spectral variability 
observed bv lAbdo et alJ d2010h for two equivalent statistical sig- 
nificance time intervals. 

The predicted light curve reaches its maximum at day 3-4, 
but the rise is not fast enough to account for the emission mea- 
sured over the first day of Fermi/LAT observations. The peak of 
the y-ray light curve actually is due to the combined action of 
the density enhancement and nova, the former providing suffi- 
cient electrons for acceleration and inverse-Compton emission 
on the photons of the latter, when the shock is still near the WD. 
As soon as the shock leaves the density enhancement and the 
vicinity of the binary, the emission drops and acceleration pro- 
ceeds in the lower-density wind. The flux decreases by a factor 
20 in the few weeks that follow the emission maximum, as ob- 
served. The predicted light curve slightly overshoots the highest 
fermi/LAT data points, and part of this deviation arises from the 
different nature of what we are comparing. The FermifLKF data 
points were actually determined using the exponentially-cutoff 
power-law fit to the time-averaged data, while our model gives 
the total lOOMeV-lOOGeV photon flux without any assump- 
tion on the spectral shape. If we derive the light curve by fitting 
the FermifLKF spectrum to our model spectra (in the 100 Me V- 
5 GeV range), we obtain ~20% lower photon fluxes because our 
model is flatter at low energies ~100MeV than what is allowed 
by the exponentially-cutoff power-law shape. Also, some data 
points are still at 2-3cr off the predicted light curve, notably at 
days 1,6, 16, and 17, and there may be several explanations for 
that: first, with a possible 1 -day offset between the outburst of the 
nova and its first optical detection, the early part of the predicted 
light curve fits rather satisfactorily the data; then, in addition to 
statistical fluctuations, these deviations may be due to inhomo- 
geneities in the density profile of the surrounding medium. The 
global no n- thermal efficiency is about ^10% after two weeks, 
and subsequently increases at a modest rate of about 1% per 
week (at least up to day 40). The X-ray thermal emission at day 
20 has a total luminosity of 5 x 10 34 ergs~' in the Swift/XRT 
band, which is a few times the measured value. When interpreted 
with the prior knowledge of the limitations of our thermal model, 
this result can actually be considered as satisfying (see Sect.|6]for 
more details). 

We emphasise that the parameters of Run 3 were not ob- 
tained through a formal multivariate fit under the constraints 
from X-ray and y-ray observations. Our aim was instead to il- 
lustrate some trends in the framework of our model, and provide 
typical values for the main parameters of the system. On that 
point, the positive influence of some sort of circumstellar den- 
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sity enhancement is noteworthy, but the statistical significance 
of the improvement brought by adding the corresponding 3 de- 
grees of freedom to the model could not be derived. The determi- 
nation of the exact shape of the CDE is beyond the capabilities 
of our model, and could more likely be done from the analysis of 
the thermal continuum and line emissi on through hydrodynam - 
ical simulations, in the same vein as lOrlando & Drake! (120121) . 
Within our model, the experimental data can be accommodated 
with a moderately more extended CDE having a lower density, or 
the opposite. The important quantity here is the swept-up mass, 
which drives the thermal X-ray output and sets the amount of 
material available for acceleration. Similarly, the orbital separa- 
tion used is still a bit smaller than the available estimates, but 
the model can accommodate slightly larger values when other 
parameters are adjusted. 

A potentially important effect in y-ray-emitting binaries is 
the anisotropic inverse-Compton scattering, and we investigated 
if that could improve the fit of our model predictions. Because 
the distribution of seed photons for the inverse-Compton scat- 
tering is not isotropic, the line-of-sight orientation with respect 
to the binary affects the observed y-ray emission. We consid- 
ered the extreme cases of superior conjunction (favouring head- 
on collisions with photons from the red giant, at least while 
some parts of the shock are still inside the binary) and infe- 
rior conjunction (favouring trailing collisions) for the system 
viewed edge-on. The results for the inverse-Compton radiation 
are shown in Fig. [9] for the two inverse-Compton components 
(the total spectra and light curves are not shown to not overload 
the plots). The impact of anisotropy on the nova component is 
modest over the first 10 days because the electron population 
is distributed nearly spherically about the source of seed pho- 
tons. At later times, the shock geometry departs from spherical 
symmetry and anisotropic effects become apparent in the light 
curve. In contrast, there is a significant impact of anisotropy on 
the RG component. The case of superior conjunction yields a 
higher flux, but the most interesting effect is on the lightcurve, 
with a flux maximum being sharper and a subsequent decline be- 
ing stronger. Inte restingly, based on an analysis of several emis- 
sion line profiles, I Shore et all d201 lb suggest that we are seeing 
the system close to superior conjunction. Yet, since the RG com- 
ponent is not the dominant one, the effect of anisotropic inverse- 
Compton on the total emission remain limited. 

6. X-ray and radio emission 

It was shown above that the y-ray emission should be interpreted 
together with the thermal emission constraints, since the latter 
dictate the geometry of the binary system and the hydrodynam- 
ics of the shock expansion. Yet, our simple model for the thermal 
X-ray emission suffers from serious limitations, which restricted 
the amount of information we could derive from the SwiftPiKT 
observations but still allowed to constrain the nova properties 
as we show below. Then, to consolidate the scenario of Run 3 
presented above, and in particular the conditions of particle ac- 
celeration, we discuss its synchrotron emission in light of radio 
observations performed in the ~l-10GHz range (see Sect. |2j. 

6.1. Thermal emission 

Our approach to compute the thermal Bremsstrahlung emission 
from shock-heated material is quite close to that adopted by 
iNelson et alJ (I2012I) . and we indeed obtained similar results to 
theirs when using the same descriptio n of the binary sys t em. A 
more sophisticated study was done by Orla ndo & Drake! d2012l) 



using multi-dimensional hydrodynamical simulations including 
radiative cooling and thermal conduction. Comparing our re- 
sults to theirs for a similar setup (for instance their best 2D 
run E44.3-NW7-CDE6.3-L40 including a CDE), it appears that 
our model for the thermal X-rays reaches the emission max- 
imum too early and largely overestimates the total luminos- 
ity over the first 2-3 weeks. At later times, our model gives 
rather good results with X-ray luminosities that are a factor ^2 
above those obtained fro m the full hydrodynamical simulations. 
Orlando & Drake! (l2012h insist on the fact that flow convergence 
on the rear side of the Mira star produces a hot and dense plasma 
that contributes predominantly to the thermal X-ray emission. 
The position of the flow convergence is what determines the rise 
of the X-ray luminosity, and this depends on the density profile 
of the ambient medium; adding a density enhancement to the 
wind density profile introduces a delay in this convergence and 
this is proposed as an explanation for the late X-ray luminosity 
maximum (Salvatore Orlando, private communication). This is 
obviously a phenomenon that our model cannot reproduce, and 
we therefore focused on the level of thermal X-rays at late times. 
In the different scenarios discussed above, we therefore verified 
that the predicted X-ray luminosity past day 20 does not exceed 
a few times 10 34 erg s _1 in the Swift/XRT band. Because of this 
shortcoming, the information about the density structure of the 
environment that is encoded in the X-ray luminosity rise is not 
accessible to us. 

If th e addition of a CDE in ou r model was motivated by the 
work o f lOrlando & Drake! d2012h . we reached different conclu- 
sions about the properties of this gas structure. The CDE of Run 
3 presented above is smaller, with a size of about 10 AU rather 
than 40 AU, and a density of 10 s cirT 3 instead of 2 x 10 6 ctrT 3 
(although our values should not be considered as precise optima, 
as already emphasised). We note that many of their assumptions 
differ from those of our Run 3, and this is probably what led 
to this discrepancy: they used an initial ejecta velocity of about 
9000 km s~', while we adopted 3000 km s~'; their ejecta mass 
is 2 x 10~ 7 Mo, which is an order of magnitude below our 2 x 
1O _6 M ; last, they considered an orbital separation of 15.5 AU 
where we used 10 AU. With such parameters, we found it almost 
impossible to account for the y-ray observations with reasonable 
parameters. Owing to the low densities, that could be only be 
done with dominantly leptonic emission. Yet, the electron ac- 
celeration efficiency would be about 10% and the electron-to- 
proton at injection close to 1. The y-ray emission maximum is 
reached early enough, but the subsequent light curve drop is far 
to o slow. In the framework o f our model, the scenario considered 
bv lOrlando & Drake! d2012l) for the thermal emission does not fit 
the non-thermal constraints. 

6.2. Synchrotron emission 

The synchrotron emission is computed for each shock element 
from the upstream magnetic field conditions and the combined 
electron population of the acceleration and cooling zones. The 
calculation includes the synchrotron self-absorption and free- 
free absorption in the shocked layer, but no absorption from 
the rest of the ambient material. In the case of Run 3, the syn- 
chrotron emission rises rapidly to reach at day 4-5 a maximum 
of about 100 mJy at 3 GHz (50mJy at 30 GHz), followed by a 
drop down to 5 mJy at 3 GHz at day 40 (1 mJy at 30 GHz). The 
time period of our synchrotron emission peak was not covered 
by radio observations, which started at day 14, but our modelled 
synchrotron emission past that date does not fit at all the radio 
data. Especially at high frequencies, the measured flux densities 
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at day 40 and beyond largely exceeds our predictions, by a factor 
of about 50 for the emission at 30 GHz. Actually, the observed 
flux densities increase with frequency, in contrast to the typical 
F(v) oc v -°- 5 of our modelled synchrotron emission. 

Both the flux mismatch and incorrect spectral distribu- 
tion suggest that the radio emission detected after the V407 
Cyg outburst is not synchrotron emission. We conjecture in- 
stead that ionisation of the red giant wind by the po werful 
UV emission during the nova outburst ( Shore et ai]|201 lb gives 
rise to thermal radio em ission with the required spectral index 
dWright & Barlowlll975l) . The time profile of the radio emission 
would then depend on the evolution of UV heating, radio absorp- 
tion, and plasma recombination. This will be addressed in more 
details elsewhere. 



7. Discussion 

In this section, we summarise the main results obtained on V407 
Cyg and discuss their implications in terms of particle acceler- 
ation as well as possible alternatives that we could not explore. 
We then consider another instance of the phenomenon, the out- 
burst of the RS Oph system in 2006. Last, we estimate the contri- 
bution of the whole population of such novae to Galactic cosmic 
rays, and the opportunities for future detections of their outbursts 
in y-ray s. 

7.1. Lessons learnt from V407 Cyg 

We aimed at modelling the Fermi/LAT observations of the V407 
Cyg outburst with the nova shock accelerating particles at a mod- 
erate efficiency compatible with the test particle approximation. 
The thermal X-ray emission measured by Swift/XKT imposes 
a relatively low density for the close environment of the bi- 
nary, which limits the proton radiation efficiency and the num- 
ber of particles that can be accelerated. In these conditions, the 
emission is mostly of leptonic origin and arises primarily from 
inverse-Compton scattering in the nova light. The observed level 
of y-ray emission can be reproduced with a non-thermal effi- 
ciency of the shock of 10% after two weeks and an electron- 
to-proton ratio at injection of 6%. The event produced protons 
(respectively electrons) with energies up to ^300 GeV (respec- 
tively ^20GeV), for a total non-thermal energy of 10 43 erg. The 
Fermi/LAT spectrum requires shock acceleration close to the 
Bohm limit in the upstream equipartition magnetic field. In or- 
der to match the observed light curve, in particular the early rise 
and rapid drop, a density enhancement on top of the stellar wind 
is required in the vicinity of the WD. Such a structure very likely 
results from accretion and orbital motion and the need for it is 
corroborated by hydrodynamical simulations and studies of the 
thermal emission time profile. We found that the light curve and 
spectrum are only mildly affected by the line-of-sight orientation 
when realistic anisotropic inverse-Compton scattering is taken 
into account. 

Overall, the V407 Cyg eruption can be understood from the 
same principles that are invoked for particle acceleration in su- 
pernova remnants (SNRs). An interesting point is that strong 
magnetic field amplification is not required, although injection 
of turbulence at the relevant spatial scales for Bohm diffusion of 
~ 1-100 GeV particles is probably needed. That difference may 
arise from the far denser upstream medium in V407 Cyg, but 
whether such a diffusion regime can effectively develop given 
the shock physical conditions would deserve a dedicated study. 

A potential issue for the scenario described above is the weak 



slowing-down of the shock and ejecta. At day 40, shock veloc- 
ities are everywhere above 2000 km s~' in our model, while line 
observations suggest that a substantial fraction of the ejected 
mass was slowed down to below 1000 km s~' at day 50. The re- 
construction of the shock kinematics from emission line signa- 
tures is not an easy task, and it may just be that we are compar- 
ing the velocities of distinct e mitting regions. Op tical spectra of 
the V407 Cyg nova reported in lShore et al.l d201 ll) show a diver- 
sity of line shapes and evolutions, which may be indicative of a 
complex underlying flow pattern with a large variety of physical 
conditions. As an example, the maximum radial velocities for the 
Ha and He II transitions differ by almost an order of magnitude 
40 days after outburst, from about 3000 to 300kms _1 for the 
blueshifted side of the lines. On top of that, part of the discrep- 
ancy may simply result from the rough prescriptions adopted for 
the shock evolution in our model, which basically neglected all 
the downstream hydrodynamics. Radiative hydrodynamics sim- 
ulations would be needed to really explore that aspect of the 
problem. Nevertheless, we investigated the effect of a more effi- 
cient shock slowing-down by forcing the velocity to drop faster 
so that it reaches below 1000 km s _1 after day 40. The most in- 
teresting consequence is that the non-thermal efficiency does not 
increase at late times and is frozen at the value it has at about 
day 10. Apart from that, the spectrum and light curve are almost 
unaffected. 

Last, we want to emphasise that the test particle approxima- 
tion led to dominantly leptonic emission and a high electron-to- 
proton ratio at injection, for the global non-thermal efficiency 
had to remain relatively low. We cannot exclude, however, the 
possibility of a comparable, if not dominant, hadronic contribu- 
tion to the emission, but that would require very efficient proton 
acceleration and our model is not suited to the study of cosmic- 
ray dominated shocks (in addition, the contribution of secondary 
electrons/positrons from charged mesons would need to be taken 
into account). For the same binary and nova parameters as in Run 
3, the FermifL AT data can be reproduced with r]i n j, P — 2 x 10~ 2 , 
r \inj,e — 2 X 10~ 4 , and f = 5, but then some parts of the blast wave 
would convert > 50% of their initial kinetic energy into acceler- 
ated particles. Although we cannot rule out that option, we note 
that the study of another nova in a symbiotic system, RS Oph, 
concluded that particle acceleration should have proceeded with 
a relatively moderate efficiency (see below). 



7.2. Another nova outburst, RS Oph 

There are many s i milarities between V 407 Cyg and RS Oph 
dAbdo et all 120101: iMunari et ail 1201 ll) . although a thorough 
analysis of emission line signatures revealed interesting dif- 
ferences betw een the V 407 Cyg 2010 event and the RS Oph 
2006 event dShore et all 1201 1). RS Oph is a symbiotic binary 
consisting of a massive WD and a RG, s eparated by about 2 AU 
and orbiting with a period of 455 days dFekel et alJ l2000). The 
system is located at a distance of 1.6kpc and was observed in 
eruption on six occasions in 1898, 1933, 1958, 1967, 1985, and 
2006. The last event in particular was covered in X-rays by 
Swift/XKT, but no y-ray instrument was available at that time to 
observe in the >100MeV range. 

We used our model to estimate the high-energy emission 
from the 2006 eruption of RS Op h. The system par ame- 
ters are similar to th o se ad opted in IWalder et al.l ([2008) and 
iTatischeff & Hernanzl d2007l) : binary separation d or b- 2AU, 
red giant radius Rrg- 150 R and mass-loss rate Mrc- 
10~ 7 M o yr _1 , ejecta mass M e j= 1O~ 6 M and initial velocity 
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V e j- 4000 km s (other parameters being similar to the base 
case for V407 Cyg). With these assumptions, the wind density 
at the RG s urface is about 10 9 cm~ 3 . The hydrodynamical sim- 
ulations by IWalder et al.l (120081) showed that the WD accretion 
and orbital motion result in a density enhancement of about 
10 cm" 3 at the WD position, with a typical scale length of 
1 AU (for comparison, we inferred for V407 Cyg a density 
enhancement at the WD of the order of the density at the RG 
surface with a scale length of 10 AU). 

In the case of particle acceleration occurring at the Bohm 
limit in the upstream equipartition field, the 2006 eruption of RS 
Oph produced protons reaching up to almost 1 TeV and elec- 
trons reaching about 30 GeV. This is quite similar to what was 
obtained for V407 Cyg in the same acceleration conditions (and 
confirms once again the robustness of the maximum particle 
energy to changes of the binary and nova parameters). Due to 
the massive accumulation of gas close to the WD, the predicted 
y-ray emission from RS Oph is dominantly of hadronic origin, 
even for electron-to-proton at injection up to 10%. For a global 
non-thermal efficiency of ^10%, the > 100 MeV flux rises up 
to 10~ 4 phcnT 2 s _1 within a day, but subsequently drops by 
more than three orders of magnitude over the following 3-4 
days. This y-ray flash actually is due to the strongly peaked 
matter distribution in RS Oph, where a relatively high density 
of target particles for pion production is found only in the 
immediate vicinity of the WD. Shock deceleration takes place 
very rapidly within the first 2 days and the level of thermal 
X-ray emission after 3 weeks is a fe w times 10 35 ergs~' in 
agreement with Sokolo ski et al.l d2006l) . The different thermal 
X-ray signatures at early times in the RS Oph 2006 event, which 
displayed powerful emission within a few days after outburst, 
and in the V407 Cyg 2010 event, which exhibited a late rise 
to a more modest emission level, very likely arise from the 
matter distribution in the innermost regions of the binary system 
(namely the density enhancement on top of the stellar wind). 

Particle accelera ti on in RS Oph was studied by 
Tatischeff & Hernanz (2007) and constrained through its 
influence on the shock structure and X-ray emission. The 
authors estimated that protons could be accelerated up to a few 
TeV in a few days but are prevented to go beyond by particle 
escape upstream of the shock. Although both approaches are 
different and we did not consider upstream escape, this is 
consistent with our result for very similar assumptions about 
the particle diffusion regime. Then, using a semianalytic model 
of nonlinear diffusive shock acceleration, the authors found that 
the post-shock temperature evolution can be accounted for if 
the shock converts up to about 60% of the instantaneous energy 
flux into non-thermal particles and is efficiently cooled through 
the escape of accelerated particles. Yet, the implied particle 
acceleration efficiency remains moderate, as in the scenario 
presented above for V407 Cyg. 

Finally, we note that there are more objects likely to be other 
examples of a fast nova explosion in a symbiotic system. The 
dramatic X-ray, optical and radio 1998 flare of CI Cam may be 
the result of a classic al nova eruption in the dense circumstellar 
medium of a Be star dFilippova et al.l l2008). More recently, the 
transient X-ray source MAXI J0158-744 detected in the Small 
Magellanic Clo ud was interpreted as a nova in a Be- WD binary 
(iLiet al.ll2012l) . 

7.3. The population of novae in symbiotic systems 

From the above elements, we can estimate the contribution of 
novae in symbiotic systems to the production of Galactic cos- 



mic rays. Their occurrence rate is in the rang e l-10yr , bu t 
not all systems are likely to be equally efficient (iLii et al.ll201 ll) . 
Assuming that the 2006 event of RS Oph and the 2010 event of 
V407 Cyg are representative, novae in symbiotic systems would 
release of the order of 10 43 erg in the form of <1 TeV protons. 
For comparison, supernova remnants are 100 times less frequent 
but are believed to release at least 10 49 erg of <1 TeV protons. 
The population of objects exemplified by V407 Cyg and RS 
Oph is therefore not expected to be a significant contributor to 
Galactic cosmic rays. 

In terms of perspective for future y-ray observations of novae 
in symbiotic systems, our work suggests that the best candidates 
would be binaries with relatively short period of order l-10yrs 
because of the high peak luminosities they can reach if they ef- 
fectively develop a significant density enhancement around the 
white dwarf (this however may be counterbalanced by a very 
rapid emission decrease for the closest systems). Optimistically 
taking our model for the 2006 RS Oph nova as a reference (be- 
cause it is brighter in y-rays), such a source could probably be 
detected by the ferm//LAT up to distances of ^lOkpc. Within 
that horizon, the occurrence rate of all novae in symbiotic sys- 
tems, irrespective of their orbital period, would be at most a 
couple of events per year (assuming a uniform distribution of 
sources and a Galactic disk of 20kpc radius). Hence, no more 
than a few events like V407 Cyg can reasonably be expected in 
the ferm//LAT remaining 5-10 years of operations. In any case, 
our analysis suggests that novae i n symbiotic syste ms are un- 
likely to be a class of TeV-emitters dAliu et al.ll2012l) . 

As we finished this work, the FermifLPX collaboration re- 
ported evidence for y-ra y emission associa ted with two classical 
novae: Nova Sco 20 1 2 dCheung et alJ2012l) and Nova Mon 20 1 2 
(Cheung et al. 2012, ATel 4310). In these cases, the companion 
is typically a low-mass main-sequence star in a tight orbit with a 
period of order of an hour. Expansion of the burning WD enve- 
lope can be accompanied by a supersonic shock as in V407 Cyg, 
so particle acceleration is also possible. The stellar radiation field 
from the companion is very weak but inverse-Compton emission 
could still occur on photons from the nova, whose contribution 
dominates the y-ray light curve in our calculations. However, 
the stellar wind is too tenuous to slow down the nova ejecta, 
which rapidly expands to engulf the whole system. Hence, lit- 
tle swept-up material can be accelerated and significant y-ray 
emission is puzzling if the environment simply consists of the 
stellar wind. One possibility is that shocks in classical novae in- 
teract with the cold, massive circumbinary discs that are thought 
to surround cataclysmic variables up to distances of several AU 
dSpruit &Taamll2001l: iDubus et al.ll2004l) . in a way that is rem- 
iniscent of the setup of Run 3 that we presented here for V407 
Cyg. A different possibility is that particle acceleration occurs 
at the temporary bow shock that forms behind the companion 
star as the ejecta expands beyond the binary orbit (but such a 
process may also be invoked in V407 Cyg and we have shown 
that it is not necessarily needed). The rate for classical novae 
in th e Galaxy is ab out 10 times higher than for symbiotic no- 
vae dLii et al.l 12006). so future ferm//LAT observations should 
be able to shed more light on this unexpected phenomenon. 
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Fig. 6. Properties of the y-ray emission for Run 1, the base 
case scenario of a shock propagating in a wind. Top panel: 
time-averaged spectrum over the first 15 days; the dashed (dot- 
dashed) curve is the inverse-Compton scattering on the red gi- 
ant (nova) photons. Middle panel: light curves in the 100 Me V- 
100 GeV band; the collective upper limit on the flux from day 
19 to 33 is 0.8 xlO~ 7 phcrrT 2 s (not shown for simplicity). 
Bottom panel: non-thermal energy in particles and cumulated 
radiation; the red and blue solid (dashed) curves correspond to 
particles (radiation); the black solid (dashed) curve is the total 
non-thermal energy (initial nova kinetic energy). 
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Fig. 7. Properties of the y-ray emission for Run 2, the scenario Fig. 8. Properties of the y-ray emission for Run 3, the scenario of 
of a shock propagating in a wind, optimized for the case of V407 a shock propagating in a wind and density enhancement around 
Cyg. Color coding and graph elements are described in Fig. [6] the WD, optimized for the case of V407 Cyg. Color coding and 

graph elements are described in Fig. [6] 
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Fig. 9. Effect of anisotropic inverse-Compton scattering. The red 
and blue sets of curves correspond to the scattering on the red 
giant and nova photons, respectively. The dot-dashed and triple- 
dot-dashed curves correspond to superior and inferior conjunc- 
tions, respectively, for the system viewed edge-on. The dashed 
curves are the isotropic case, shown for reference. Injection frac- 
tions and diffusion efficiencies are the same as in Run 3. 



18 



